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This letter contains several formulas and algorithms involving the (affine) Weyl groups, 
which play central roles in the author's present work. Though some are surely known to 
other experts, some are new, and they are all collected (and analysed) here to benefit also 
i) ". the researchers and students less experienced with these matters. An indication of the 



^ ' need for good algorithms is provided by the number of papers (see e.g. [1-4]) on it. 

One place the finite Weyl groups arise is in the modular transformation [5] 

K 

ci\ xx{-'i-/r) = ^ Sxt^Xt^{r) (la) 

of affine characters, where the sum is over the level k integrable highest weights fx G 

for the given affine Kac-Moody algebra Xr^\ for the k satisfying A G The level /c is a 
nonnegative integer, and the finite set is defined in terms of the colabels by 

= {(Ai, . . . , A,) I a^Xi + ■■■ + a^Xr < A;, A, G Z>o} (16) 

We identify the r-tuple with the combination AiAi -(-■ ■ ■ + XrAr in terms of the fundamental 
weights Aj. An expression for this unitary and symmetric matrix is [5] 

S,, = V (det«;) exp[-27ri + P) " (/^ + ] (i,) 



1 



where p = (1, . . . , 1) and k = k + . The colabels o^, the finite Weyl group W, the 
number s = ill^+ll / ^||(5^*/Q^||, the dot product, and the dual Coxeter number /i^, wiU 
be exphcitly given below for each algebra. 

The challenge raised by several researchers is to find an effective way of computing 
this matrix S. For instance, in modular invariant classifications [6], literally millions of 
S entries typically must be computed for each algebra. The difficulty is that the Weyl 
group W — though finite — can be very large and unwieldy. This is particularly acute 
for the exceptional Xj,: the Weyl groups of G2, Eq, Ef, Eg respectively have 12, 1 152, 
51840, 2 903 040, and 696 729 600 elements. We wiU show that this sum (Ic) is easy to 
evaluate for the classical Xr^^ , and so a way to evaluate it for the exceptional ones would 
be to find a large classical Weyl subgroup. Our algorithm involves computing 0((log|VF|)^) 
trignometric functions and performing 0((log|VF|)^) arithmetic operations. 

One of the many places in which the matrix S appears is Verlinde's formula [7] for 
the fusion coefficients N^^. Write for the weight (0, 0, . . . , 0) e P^. Then 



Verlinde's formula is not of much direct computational value: to leading power in k, it has 

k"- 

il ■ ■ ■ Uj. 



terms. This also gives the size of S. More effective is the Kac- Walton formula [8,9], 
which expresses the fusion coefficients in terms of X^ tensor product coefficients T^^ = 
multi,(LA ® L^): 

^Am= E(d^t^)^A/ (2c) 

where W is the affine Weyl group; its action on the weights Pj; will be discussed in §2.2. 

Background material on simple Lie and affine algebras can be found in [8,10-12]. For 
an account of the rich role of the matrix S and fusion rules in rational conformal field 
theory and elsewhere, see e.g. [13,14] and references therein. 

In this paper we focus on three tasks: effective formulas for computing 5" entries; fast 
algorithms for finding the preferred affine Weyl orbit representative [A] — w.X & P^; and 
super-fast formulas for obtaining the parity e(A) = det('u;) for that representative, when 
the weight w.X itself is not needed. In Section 4 we sketch some applications. 

2. Review and statement of the problems 

Introduce the notation Pj, for all XiAi = (Ai, . . . , Ar-) E 7/". It inherits the usual dot 
product of P^. When the level k is understood, we'll write 
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Also, we'll abbreviate 'algebra Xr^^ at level A;' by 'Xj.,^'. We will write '[yj' for the largest 
integer not greater than y — e.g. [ttJ = 3 = [3J. As always, k = k + hy . 

2.1. Nice weights and useful symmetries. If either A or is a 'nice' weight, then significant 
simplifications to (Ic) can be made. In this letter we are interested rather in generic 
weights A, ^ in (Ic) and so we will not devote much space reviewing these. 

The denominator identity of Lie algebras permits us to express the sum (Ic), whenever 
A is a multiple of p, as a product over the positive roots: 

n -r/2\ I TT <^ • / ("^+l)(At + P)-«N /o\ 

Smp,^. = Ki '^'^\s\ II 2sm(7r^ ) (3) 

aeA+ 

The most important choice is m = 0, which yields a product formula for the quantum- 
dimensions Sxq/Sqo- For a rank r Lie algebra Xr, the Weyl group has approximate size 
r! ^ y/27rr{r/eY, while the number of positive roots grows like r^, so (3) is a significant 
simplification. 

Remarkably, we will find that all entries of S are about as accessible as Sxo, in par- 
ticular, via equations with complexity growing like r^. 
An important observation from (Ic) is that the ratio 

|^ = chA(exp[-27ri^]) (4a) 

is a finite Lie group character evaluated at an element of finite order [3] . Effective expres- 
sions exist for some ch^, implying that the corresponding entries Sx^ can be computed 
quickly. The weight multiplicities mult2:(A)(/5) for simple algebras of small rank and small 
weights are given explicitly in [15]. When L(A) has few dominant weights and each 
of these has a simple Weyl orbit, then (4a) can be evaluated explicitly. For instance, for 
Xr = A^, when A is a fundamental weight A^, we get 

^^|^=exp[27rimt(p + p)/(r + l)Av] ^ Xi,---Xi^ (46) 

where t(z/) = Yli=i^^i Xi = exp[— 27rip"'"[i]/K;] for orthogonal coordinates i'~^[i] = 

Another application of (4a) is for the use of branching rules in replacing Sx^/Sq^ for 

Xr^^ with a sum (over A') of S'y^,/SQ,^, for a simpler algebra Xi.'^'^^ of equal rank. E.g. the 
embedding C Eg, has the branching rule (A2) -|- (Ag) = (Ai) (see e.g. [30]), and hence 
the ratio Sa-^^/Sq^ equals the sum S'j^, ^, / S'q,^, + S'j^, for an appropriate choice of 

p' (depending only on p). Later we will work this out explicitly, for both D4 C F4 and 
A2 C G2. 

Related to this, for fixed 7 G P+, the ratios (4a) form a one-dimensional representation 
of the fusion ring. That is, 

•^^7 '^M7 _ Art/ ^1^1 /r^ 
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Using (2c), we can obtain S entries for 'larger' weights from those of smaller weights. 

Considerable simplifications [3] to the calculation of chA(exp[— 27ri ^]) occur when 
some Dynkin labels of fi are 0. While this is not of direct value to us (because of the 

in (4a)), it can be very useful as an approximation when the level k is large (see e.g. 



Incidentally, the matrix 5" obeys several symmetries. Best known are those related 
to the symmetries of the extended Dynkin diagrams (see e.g. [12]). Thus determining 
one S entry automatically yields several others. For example, consider = A^. De- 
fine an involution C by C(Ai, . . . , Aj.) = (A^., . . . , A2, Ai), and an order-(r + 1) map J by 
J(Ai, . . . , A^) = (Ao, Ai, . . . , Xr-i) where Xq = k — Yll^i Aj. Then both C and J are per- 
mutations on P^, corresponding to the symmetries of the extended Dynkin diagram of A^. 

The various powers J* are the simple- currents of A^^K We have Scx,i^ — '^x a 



Similar equations hold for the other algebras (see e.g. [12]). 

Another symmetry — more complicated but more powerful — is the Galois symmetry 
[16] discussed in §4.4 below. Making it more accessible is a big motivation for the paper. 

For the classical algebras, another symmetry (rank-level duality [17]) of the S matrix 

tells us that small level acts like small rank. In particular, the S matrix for sl(n)^ and 

sl(A;)^ are closely related, as are so(n)^ and so(A;)^, and Cr,k and Ck,r- 

Finally, simplifications (called fixed-point factorisation [18]) occur whenever A or // is 
a fixed-point of a (nontrivial) simple-current — the simple-currents J* for A^^ are given 
above. If 99 is fixed by J^, then for any weight of A of sl{n)k we can write Sip\ as a product 
of n/d ^'-entries for sl{d)kd/n- For a very simple example, take sl{n)k when n divides k: 
the unique J-fixed-point is = (^, . . . , ^), and Sx^^p — ±(^)*^"^~^)/^ or 0, depending on A. 

Incidentally, an analogue of fixed-point factorisation works for the symmetries (called 
conjugations) of the unextended Dynkin diagrams. For example, consider the Zs-symmetry 
triality of D^: a weight /i invariant under it obeys jJ-i = fis = jJ-^, and the D^^ ratio Sxh/Sqh 
for such a /i equals a sum '^A'ai'/'S'o>' of G2^^ ratios, where ®(A') = (A) are G2 C -D4 
branching rules (see e.g. [30]), and ^' = (/X2,3/xi + 2) e Pf+^(G'2). Applications of these 
fixed-point factorisations is given in [31]. 

In practice the remarks of this subsection are quite effective. For example, consider 
^3,6- There are 7056 S entries: any involving J-?0, J-?(l, 1, 1), J-?(2, 2, 2), J^A^, J-?(l, 0, 1), 
J-'(l, 1, 0), (2, 0, 0), or any J^-fixed point, are immediately calculated from the above. 
Finally, exploiting the Galois symmetry (see §4.4 below), we reduce the calculation of the 
7056 entries to precisely 1 less pleasant calculation: that of -S^a 0,0), (3,0,0)- Of course, 
for large k and r, most weights won't be 'nice' in this sense, and another approach is 
required. 

2.2. Calculating Weyl orbits. The affine Weyl group W is generated by the refiections rj 
through simple roots: explicitly, for any A G Pz, the Dynkin labels of r^A are 



[6]). 



SJx,^^ = exp[27ri t(/i)/(r + 1)] Sx^, 



(6) 



iji^j — Xj — XiAji 



(7a) 
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where Aji are entries of the Cartan matrix of Xr . It is important that W can be expressed 
as the semidirect product of the finite Weyl group with the coroot lattice interpreted as 
an additive group — i.e. each w can be uniquely identified with a pair (w, a), for w E W 
and a e , and w\ — w{X) + na. It is often more convenient to express this action using 
the notation 

w.X:=w(X + p)-p {7b) 

Given any level k weight A = (Ai, . . . , A^) e Pz, there are two possibilities: either A + p is 
fixed by some affine Weyl reflection w, i.e. w.X = A, and we shall call A null; or the orbit 
W.X intersects the fundamental alcove in precisely one point: 

[A] := W.X e (7c) 

Define the parity e(A) to be or (iet{w), depending on whether or not A is null. This 
preferred orbit representative [A], and/or the parity e(A), are often desired — see §4. 

A very useful fact is that A will be null if A^ = — 1 for any i = 0, 1, . . . , r. 

A method for finding [A] and e(A) is proposed in [12]. Put e — +1 and jj, = X. 

(1) If > for all z = 0, 1, . . . , r, then [A] = ^ and e(A) = e. 

(2) Otherwise, let < i < r be the smallest index for which /x^ < 0. Replace /x with r^./x, 
and e with — e. Goto (1). 

It does not appear to be known yet whether their method will always terminate, 
although in practice it seems to. We give an alternative next section. 



3. The formulas and algorithms 

The essence of the following S matrix formulas is the observation that an alternating sum 
over the symmetric group also occurs in determinants. The point is that determinants can 
be evaluated very effectively, using Gaussian elimination, and is easy to implement on a 
computer. By comparison, determinants are much more accessible than permanents. This 
suggests that we look no further than (Ic) for an effective algorithm. 

The classical 5' matrices all reduce to evaluating one or two determinants. Our strategy 
for the exceptional algebras is to find a classical Weyl group of small index. For example 
W{D^) is contained in W{E^) with index 135. Once we find coset representatives for the 
135 cosets in W{D%)\W{Es), then the Es S matrix would be the alternating sum of 135 
Dg S matrix entries. Equally important, these coset representatives also permit effective 
algorithms to find [A] and e(A), as we will see. 

A priori, it could be hard to find these coset representatives, but here the task is made 
elementary by the following result [19] (there's a typo in their definition of -D*): 

Theorem. Let ^ he a subsystem of a root system $. Choose simple roots a & J 
for \1/ which are positive roots $+ in $ (this is always possible). Define = {w G 
II^($) \ w{J) C Then every element o/II^($) can be uniquely expressed in the form 

dw' where d e and w' e W{^). Furthermore, d is the unique element of minimal 
length in the coset dW{^). 
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In particular, the set of all d~^, for d G -Dijf, constitute the desired coset representa- 
tives for W{'if)\W{^). We chose the largest classical subsystem possible for each of the 
exceptional root systems. The results are explained in the following subsections. 

To give the extreme example, we read off from Table 2 that the W {Ds)\W [Eg) coset 
representative we call C135 is the following composition of simple reflections: 



r7r6r5r4r3r2rir8r5r4r3r2r6r5r4r3r8r5r4r6r5r8. 



It acts on a weight A G PiEg) as follows: write u = C135.A, then using notation defined in 
§3.E8 we get 



1 
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Root system data is explicitly given in e.g. [10,12]. For convenience we reproduce it 
in the following subsections. 

The best way to determine the parity e(A) of §2.2 seems to be to use (3): 

e(A) = n sign{sin(7r iA±^)} (8) 

The (easy) proof that (8) holds is that the product of sine's can be expressed as alternating 
sums over the Weyl group, thanks to equation (3), and that when evaluated at any weight 
A G they yield a positive value (proportional to the quantum-dimension). Equation (8) 
has been kicking around for years — see e.g. the authors earlier work, and also [20]. Later 
in this section we make significant improvements to (8) for all algebras except Ai^\ 

We often require ordering lists of numbers (lists whose length is approximately 
the rank r of the algebra). Straightforward ordering algorithms have a computing time 
growing like A^^, but more sophisticated algorithms have order log A^, and so are relevant 
for large rank. See [21] for details. 

Many of our formulas involve taking determinants of N x N matrices, where A^ ap- 
proximately equals the rank. Using Gaussian elimination, this involves approximately A^^ 
operations. We will briefly discuss the error analysis in §4.1 — see also [22]. 

Some of the following (namely the S matrix for the classical algebras) has appeared 
in [17], and B. Schellekens has used but not published similar S matrix formulas (see his 
webpage at [http: / /www. nikhef.nl/ ~t58/kac.htm]|) . Determinant formulas for the charac- 
ters of the classical Lie algebras go back to Weyl and are discussed in §24.2 of [23]. Other 
approaches to computing generic S entries can be found in [1,2,4]. 
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3. A. The A-series. The colabels for Xr = Ar = sl(r + 1) all equal 1. The scale factor in 
(lb) is s = 'f('^+'^y^/y/r + 1, and k = k + r + 1. To the weight A = (Ai, A2, . . . , A^) associate 
the orthogonal coordinates X[£] = J2l=e for £ = 1, . . . , r + 1, and write t(A) = J^g^i 
and A"'"[j] := (A + p)[j] — r + 1 — j + A[j]. The finite Weyl group is the symmetric group 
Sr+i, acting on A by permuting the components X[£]. The affine Weyl group also acts on 
a level k weight A by translation: A 1-^ A + A;q: for any a belonging to the Ar root lattice, 
i.e. a= {ci,..., Cr+i) G Z''+^ Cj = 0. Then 

K [V ~T' i- I 

'det' in (A.l) denotes the determinant of the (r + 1) x (r + 1) matrix whose (i, j)-th entry 
is provided. 

The fastest algorithm known to this author for finding the orbit representative [A] of 
(7c) is as follows. 

(i) Write Xi = X+[£]=r + l-e + Aj for £ = 1, 2, . . . , r and put Xr+i = 0. 

(ii) Reorder these so that xi > X2 > ■ ■ • > Xr > x^+i > 0. If xi — x^+i < k, then done; 
goto step (iv). 

(iii) Otherwise for each i = 1, . . . , (r + l)/2 put mj = [(xj — Xr+2-i) / K'I', if is even 
let Xi = Xi — miK/2 and Xr+2-i = Xr+2-i + miK/2, while if rrii is odd let Xi = 
Xr+2-i — {j^i + 1)k/2 and = Xi + {rrii + 1)k/2. Return to step (ii). 

(iv) The Dynkin labels [A]j of the desired weight are 

[A]i = yi- Vi+i - 1 

The weight A is null iff this algorithm breaks down: i.e. two (or more) x^'s are equal. 
The parity e(A) is the product of the signs of the permutations in step (ii), together with 
(— l)"^'s from step (iii). An alternative is provided in [24]. 

Equation (8) becomes 

6(A) = sign TT sin(7r^^H^^) (A2) 

'i-<i<j<r+l 

3.B. The B-series. For = — so(2r + 1), a( — — 1 and all other colabels equal 2. 

2 

The scale factor is s = i*^ /2, and n = k+2r—l. To the weight A = (Ai, A2, . . . , A^-) associate 
the orthogonal coordinates X[£] = Yl^Zt + Ar/2 for £ = 1, . . . , r (put J2lZr — 0) 
write X'^[j] := (A + p)[i] — r + ^ — j + X[j]. The finite Weyl group is a group of order 
2^ ■ r!, acting on A by arbitrarily permuting and changing signs of the components X[i]. 
The affine Weyl group also acts on a level k weight A by translation: X ^ X + ka for any 
a belonging to the Dr root lattice, i.e. a = (ci, . . . , c^) G with Ci is even. Then 

Sx^ = «--/22-ii-'--det(sin(27r^^il^))i<,,,<, (S.l) 
The best algorithm known to this author for finding the orbit representative [A] is: 
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(i) Write xi = \+[e\ = r-i+\ + + for £ = 1, 2, . . . , r. 

(ii) By adding the appropriate multiples of k to each xg, find the unique numbers yi, . . . ,yr 
such that X£ = yg (mod k) and —k/2 < ye < k/2. Let m be the total number of k 
added: i.e. if Xi — yi + rriiH, then m — Yll^i rrii. We need to know later whether m is 
even or odd. 

(iii) Replace each yi with its absolute value. Reorder these so that j/i > y2 > • • • > J/r > 0. 

(iv) If m is odd, replace yi with k — yi. 

(v) Then the Dynkin labels [A]j of the desired weight are 

=yi- yi+i - 1, for 1 < z < r, and [Aj^ = 2yr - 1 

The weight A is null iff two of the t/j's in (iii) are equal, or if at least one equals 0. 
The parity e(A) is (— l)"* times the sign of the product Ylyi in (ii), times the sign of the 
permutation in (iii). 

Equation (8) becomes (note the decoupling of the ith and jth terms here) 

e(A) = sign n sin(7r^^) ff {cos(27r ^^) - cos(27r ^^^)} (S.2) 

l<i<r l<i<j<r 

where we used the identity 2 sin(a; — y) sin(a; + y) = cos(2j/) — cos(2a;). 

3.C. The C-series. The colabels for Xr = Cr = sp(2r) all equal 1. The scale factor is 
s = 1^ /2''/^, and k = k + r + 1. To the weight A = (Ai, A2, . . . , A^) associate the orthogonal 
coordinates X[i] = J2l=i A, for £ = 1, . . . , r, and write X+[i] = {X + p)[i] = r + 1 - i + X[i] . 
The finite Weyl group is a group of order 2^ • r!, acting on A by arbitrarily permuting 
and changing signs of the components X[£]. The affine Weyl group also acts on a level 
k weight A by translation: A 1— > A + A;q; for any a belonging to (2Z) ® • • • ® (2Z), i.e. 
a = (ci, . . . , Cr) e Z*", where each Cj is even. Then 

= (2/«)^/2i'"'-Met(sin(7r ^iMi^))i<,,,.<, (C.l) 

The best algorithm known to this author for finding the orbit representative [A] is: 

(i) Write xe = X+[i]=r + l-i + 

(ii) By adding the appropriate multiples of 2k to each xg, find the unique numbers 
yi, . . . ,yr such that X£ = yg (mod 2k) and —k < ye. < k. 

(iii) Replace each j/j with its absolute value. Reorder these so that k > j/i > j/2 > • • • > 

yr > 0. 

(iv) Then the Dynkin labels [X]i of the desired weight are (put j/^+i = 0) 

[X]^ = yt- yi+i - 1 

The weight A is null iff some yi = 0, or two j/j's are equal. The parity e(A) is the sign 
of the product YlUi (^0; times the sign of the permutation in (iii). 
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Equation (8) becomes (note the decoupling of the ith and jth terms here) 

e(A)=sign n sin(7r^^) FT {cos(7r - cos(7r ^^)} (C.2) 

l<i<r l<i<J<'" 

3.D. The D-series. For = D^- = so(2r), = a^-i = = ^ all other colabels 
are 2. The scale factor in (lb) is s = i''('"~i)/2 and k = A; + 2r — 2. To the weight 
A = (Ai, A2, . . . , Ar) associate the orthogonal coordinates X[£] = X]i=£ -^i ^ — '"~2 '^"^ 
£ = 1, . . . , r {J^lZl = 0), and define A+[£] = (A + p)[£] = r-£ + A[£]. The finite Weyl group 
has order rl2^~^, acting on A by permuting the components \[£] and changing an even 
number of signs. The action of the affine Weyl group includes the translations A 1— ^ A + /cct 
for any a belonging to the Dj. root lattice, i.e. a = (ci, . . . , Cj.) e Z*^ where Cj is even. 
Then 



Sx, = ^-r/^2^-^irir-i) |det(cos(27r ^J^l^))i<,,,<. (D.l) 

+ (-i)'-det(sin(27r ^^ili^))i<,,,.<,} 



The best algorithm known to this author for finding the orbit representative [A] is: 

(i) Write xe = X+[i]=r-i + YJiZl K + ^^4^. 

(ii) By adding the appropriate multiples of k to each xi, find the unique numbers yi, . . . , 
such that Xi = yi (mod k) and —k/2 < ye < Let m be the total number of k 
added: i.e. if Xi = yi + rriiK, then m = ^[^^ rrii. We need to know later whether m is 
even or odd. 

(iii) Let m' = ±1 be the sign of Yii Hi (take m' = 1 if the product is 0). Replace each y^ 
with its absolute value \yi\. Reorder these so that k/2 > yi > y2 > ■ ■ ■ > yr > 0- 

(iv) Replace yr with m'yr — ^yr- If w is odd, replace yi with k, — yi and yr with —yr- 

(v) Then the Dynkin labels [A]j of the desired weight are 

[A]i = - yj+i - 1 l<i<r-l, [X]r = yr + yr-i - 

The weight is null iff two j/i's in (iii) are equal. The parity e(A) is the sign of the 
permutation in (iii). 

Equation (8) becomes (note the decoupling of the ith and jth terms here) 

e(A)=sign n {cos(27r ^iiS) - cos(27r ^^)} {D.2) 

l<i<j<r 

Curiously, this sign has the interpretation as the sign of the permutation which orders the 
numbers cos(27r A"'"[z]/k) in increasing order. 
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3.E6. The algebra Eq. The colabels for Eq are 1,2,3,2,1,2, respectively. The scale factor 
in (lb) is s = l/\/3, and k — k + 12. To the weight A = (Ai, A2, . . . , Ae) associate the 
orthogonal coordinates (interpret ^^^=2 ^ ^) 

m = + ^ A, for £ < 5, A[6] = ^h±^ + 3A3 + 2A4 + A5 + 2Ai 

The 27 coset representatives for W{D5)\W{Eq) are given in Table 1, recursively de- 
fined in terms of the reflections r^, z = 1, . . . , 6, through the simple roots ctj. Then 



Sx^ = f;(det ce) { det(cos(27r ^e6.1) 

- idet(sin(27r ^ ' ^ ))i<r,j<5} 

Recall that the numbers (c£.A)"''[z] are the orthogonal coordinates of the weight q(A + p). 
The orbit representative [A] in (7c) can be found as follows. 

(i) Choose one of the 27 coset representatives q. Write Zi := (q.A)"'"[6 — i], z = 1, . . . , 5, 
and zq = (c£.A) + [6], for the orthogonal coordinates (reordered) of q(A + p). (The 
'[6 — i]' is chosen to make the first five indices here consistent with those of -D5.) 

(ii) Add an appropriate multiple 3nK/2 of 3k/2 to zq to get < xq < 3k/2; replace 
xi = nK/2 + zi, and Xi = —nK/2 + for 2 < z < 5. Use the D5 algorithm for the 
present value of k to get jj/s] < < Us < y2 < Vi and yi + y2 < 

(in) Write t = {-yi - ^2 - 1/3 - 1/4 + 1/5 + l/6)/2. If t > and yi + 1/2 H h ye < 2k, then 

the Dynkin labels of the desired weight [A] are 

[A]i = t-l , [X]e = 1/4 + ^5-1 , [A]i = ye-i - Vr-i - 1 2 < z < 5 

(iv) Otherwise try the next q in (i). 
Equation (8) becomes 

e(A) = sign n {cos(27r^^) - cos(27r^^)} {E6.2) 

l<i<i<5 

yj, . siA+[l]+A+[5], , A+[6]+S2A+[2]+S3A+[3]+s'A+[4] 
X {cos(7r — !^)-cos(7r — — — — —)} 

Si 

where the second product is over all 8 choices of signs si, S2, S3 G {±1}, and s' = S1S2S3. 

3.E7. The algebra Ef. The colabels for E-j- are 1,2,3,4,3,2,2, respectively. The scale factor 
in (lb) is s = — i/-s/2, and n = k + 18. To the weight A = (Ai, A2, . . . , Ay) associate the 
orthogonal coordinates X[£] = Aiz^ _)_ Y^j^^ Aj for £ < 6, (interpret Xli=2 -^i ^ ^) 
A[7] = 2Ai + 3A2 + 4A3 + 3A4 + 2A5 + Ae + 2A7. 
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The 63 coset representatives for W{AiDq)\W{E7) are given in Table 1, recursively- 
defined in terms of the reflections r,, z = 1, . . . , 7, through the simple roots ctj. Then 



Sxf^ = y](det a) sm(7r ^ ^ ^ M L J ^ ^^^-^^ 



1=1 

X {det(cos(27r ^ ^ ; ^J^ )^^..^^ _ det(sin(27r ^ ^ ; ^ i<iJ<6)} 

hi hi 

Recall that the numbers (Q.A)+[i] are the orthogonal coordinates of the weight q(A + p). 
The orbit representative [A] in (7c) can be found as follows. 

(i) Choose one of the 63 coset representatives q. Write Xi := (q.A)"'"[7 — z], z = 1, . . . , 6, 
and xr := (q.A)+[7], for the orthogonal coordinates of C£(A + p). (The '[7 — z]' is to 
make the indices here consistent with those of Dq). 

(ii) By adding an appropriate multiple of 2k to and replacing if necessary by k — x^, 
get < j/y < 3k/2. Similarly, using the Dq movements for the given value of get 
< lye] < 2/5 < 2/4 < 2/3 < y2 < yi and yi + y2 < «• 

(iii) Write t = {—yi — 2/2 — 2/3 — 2/4 — 2/5 + 2/6 + 2/7)/2- If both t > and t e Z, then the 
Dynkin labels of the desired weight [A] are 

[A]i = t-l , [Ajy = 7/5 + ^6-1 , [A]j = X7-i - X8-i - 1 2 < Z < 6 

(iv) Otherwise replace yr = k - y7, yi ^ y^ - k/2, ye ^ yi + k/2, and yi = k/2- yr-i for 
z = 2, 3, 4, 5. Define t as in (iii); if both t > and t E Z then write [A]j as in (iii). 

(v) Otherwise try the next q in (i). 

Equation (8) becomes 

e(A) = sign sin(7r^i£l) ff {cos(27r^^) - cos(27r^^)} {E7.2) 

l<i<j<6 

^rts ( giA+[l]+g2A+[2] + a3A+[3]+A+[4] ^ A+[7] + asA+fS] - 8^A+[6] 
X {cos(7r — — — —) — cos(7r — — — —)s 

Si 

where the second product is over all 16 choices of signs si, 52,^3,55 e {±1}) and s' = 

S1S2S3S5- 

3.E8. The algebra Eg. The colabels here are 2,3,4,5,6,4,2,3, respectively. The scale factor 
in (lb) is s = 1, and k = k + 30. To the weight A = (Ai, A2, . . . , Ag) associate the 
orthogonal coordinates X[£] — ^^s^^ + Ylt=8-e^i for 1 < £ < 7, (interpret ^^=7 Ai as 0) 
and A[8] = Ai + 2A2 + 3A3 + 4A4 + 5A5 + fAe + 2A7 + fAg. 

The 135 coset representatives for W{Ds)\W{Es) are given in Table 2, recursively 
defined in terms of the refiections r^, z = 1, . . . , 8, through the simple roots ctj. Then 

^A,. = ^ f:(det c) {det(cos(27r ^'''^^^^^ ^^^^^ ))i<,,<8 {E8.1) 

+ det(sin(2. "='-^'"H^'"'^'i )).<,,<,} 
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Recall that the numbers {ci.X)~^[i] are the orthogonal coordinates of the weight q(A + p). 
The orbit representative [A] in (7c) can be found as follows. 

(i) Choose one of the 135 coset representatives q. Write Xi :— (q.A)"'"[9 — i], i = 1, . . . , 8, 
for the orthogonal coordinates of q(A + p) (the '[9 — i]' is to make the indices here 
consistent with those of Ds). 

(ii) Using the Dg movements for the given value of k,, get < Ivsl < 1/7 < ■ ■ ■ < Hi and 
yi + y2< «• 

(iii) Write t = yi — 1/2 — 1/3 — y4 — 1/5 — ye — 1/7 + Us- If ^ > then the Dynkin labels of the 
desired weight [A] are 

[A]7 ^t/2-1, [\]s = y7 + y8-l, [A]^ = y^+l - 2/^+2 - 1 1 < i < 6 

(iv) Otherwise replace each j/j with k/2 — yg-i- Compute the resulting t as in (iii); if t is 

positive than write [Aj^ as in (iii). 

(v) Otherwise try the next coset representative q in (i). 

Equation (8) becomes 

= sign n {cos(27r^ifel) - cos(27r^^)} {ES.2) 

l<i<j<8 

yj. . siA+[l] + S2A+[2]+53A+[3]+A+[7], , S4A+K] + S5A+[5] + s'A+[6] + A+[8] , 
X {cos(7r — — — —) — cos(7r — — — —) 

Si 

where the second product is over all 32 choices of signs si, S2, S3, S4, S5 e {±1}, and where 

s' = S1S2S3S4S5. 

3.F4. The algebra F4. The colabels are 2,3,2,1, resp. The scale factor in (lb) is s = 1/2 
and K, = k + 9. To the weight A = (Ai, A2, A3, A4) associate the orthogonal coordinates 

(A[l], A[2], A[3], A[4]) = (Ai + 2A2 + ^Xs + A4, Ai + A2 + ^Ag, A2 + ^Ag, ^Ag) 
Let Tj, i = 1, . . . , 4 be the reflections through the simple roots cuj. Define the matrix 



1 

'=2 



(\ 1 1 -1^ 
11-11 
1-11 1 
\l -1 -1 -1/ 



It is orthogonal, with determinant +1. Then 

8 rn f . fc. . f . fc. (c-A)+[i]//+[j],, 
Sx,i. = ^{det(sin(27r^^^^^^))i<ij<4 + det(sin(27r^ ))i<m<4 

+ det(sin(27r ^"'-^^^['J^^[^'^ ))i<,,,-<4} (F4.1) 
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In (F4.1), (c.A)"''[z] denotes the ith. orthogonal coordinate of the matrix product of c with 
the column vector with entries (A + p)[j], and c* means the transpose(=inverse) of c. 
The orbit representative [A] in (7c) can be found as follows. 

(i) Let Xi = X'^[i] for i = 1,2,3,4. Apply the B4 algorithm (steps (ii)-(iv)) with the 
given value of k to obtain the 4-tuple {i/i), n > yi > 1/2 > Us > Va > 0- 

(ii) If yi > y2 + Vz + Va, proceed directly to (iii); otherwise first replace yi with the 
coordinates of the matrix product of c* with the column vector (j/j), and if necessary 
change the sign of the resulting 1/4 so that it's positive. 

(iii) The desired weight [A] has Dynkin labels 

[A]i = y2 - 2/3 - 1 , [A]2 = - 2/4 - 1 , [A]3 = 2^4 - 1 , [A]4 = yi - 2/2 - 2/3 - 2/4 - 1 
Equation (8) becomes 

6(A) = sign nsin(7r^) J] {cos(27r^) - cos(27r^)} 

i=l l<i<j<4 

n, , A+[4], , A+[l] +S2A+[2] +S3A+[3],, , 
{cos(7r — ^) - cos(7r — ^-^ ^)} (F4.2) 

S2,S3 

where the product is over the four choices of signs S2, S3 G {±1}- 

The embedding S4 C F4 (resp. C F4) lets us write the ^4^^ quantities in terms of 
the more familiar S4^fc_|_2 (resp. -D4,fe+3) ones, and this can be useful (see [31]). In particular, 
a weight A e P\{F/^ corresponds to the orbit of A' = (A2 + A3 + A4 + 2,Ai,A2, A2 + 
A3 + 1) e P^^ip^^ under P' : A3 <-> A4 (which changes the chirality of D4 weights) 
and triality T' : Ai — > A3 — > A4 — > Ai (we prime the D4 quantities). We get S\^J, = 

y i=o j=o 

where by (7') = (A) are D4 C F4 branching rules (see e.g. [30]). 

3.G2. The algebra G2- The colabels here are 2,1, respectively. The scale factor in (lb) is 
s = — 1/a/3, and n = k+A. To the weight A = (Ai, A2) associate the orthogonal coordinates 
(A[l], A[2], A[3]) = (Ai + A2, Ai, — 2Ai — A2). The finite Weyl group is the symmetric group 
^3, together with the sign change A 1-^ —A. Then 

Sx,ij. — I- {€{2771771' +mn' +nm' + 2nn') + c{—mm' — 2mn' —nn' +nm') (G2.1) 
v3k 

+ c{—mm' + mn' — 2nm' — nn) — c{—mm' — 2mn' — 2nm' — nn') 
— c{2mm' + mn' + nm' — nn') — c{—mm' + mn' + nm' + 2nn')} 

where for convenience we put c{x) = cos(27ra;/3K), m = A"'"[l], n = A"''[2], m' = IJ''^[1], and 
n' = fi+[2]. 
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The orbit representative [A] in (7c) can be found as follows. 

(i) Let Xi = A"*" [i] . 

(ii) By adding the appropriate multiples of k to each Xi, find the unique numbers < 
yi, ^2, < K such that Xi = Ui (mod k). 

(iii) If 2/1 + 2/2 + 2/3 = 2k, then replace each 2/i with n — yi. 

(iv) Reorder the 2/i so that 2/3 > J/i > 2/2- The Dynkin labels of the desired weight are 
[A]i = 2/2-1 and [A]2 = 2/1-^2-1- 

Equation (8) becomes 

e(A) = sign{s(3a) s{h) s(3a + h) s{2,a + 26) s{a + 6) s(6a + 36)} {G2.2) 

where for convenience we write a = Ai + 1 and 6 = A2 + 1, and put s{x) = sin(7ra;/3K). 

The embedding A2 C G2 lets us write the 6*2, fc quantities in terms of the more familiar 
^2,fc+i ones, and this can be useful (see [31]). In particular, any weight A G P^(G2) 
corresponds to a pair A' = (Ai,Ai+A2 + l),C"A' = (Ai+A2 + l,Ai) of weights inPj+i(A2). 
We get Sxiu. = i S'^^, - i S'^^, and 

7' 

where 6y (7') = (A) are A2 C G2 branching rules (see e.g. [30]), and C : iJ,[ ^ ji'2 is 
A2 charge-conjugation. 

4. Further remarks 

Error analysis and exact results. The error analysis of the Gaussian elimination 
method for computing determinants is surprisingly subtle and has been the subject of 
extensive study — see e.g. [22] . It depends on the size of the pivots, but very typically the 
error grows by a factor on the order of -Jn for an n x n matrix. A good idea however is to 
make sure in your program that the pivots are never too close to 0. 

Because our quantities Sx^ are all cyclotomic integers (up to a global rescaling) in 
some field Q[exp[27ri/n]], the obvious way to make all calculations exact is to do it over the 
polynomials p{x) with integer coefficients. The desired numerical value would then simply 
be p(exp[27ri/?i]). These polynomials add, multiply etc in the usual way, but we regard 
two polynomials as equal if they differ by a multiple of the cyclotomic polynomial (f)n{x). 
The cyclotomic polynomial is the polynomial of smallest positive degree, with integer 

coefficients, which has exp[27ri/n] as a root. For example, </>p(a;) = a;^~^ + a;^~^ + hi 

and (f)nix) = 0p^...p( (a;"/^^ '^*) where n = ■ ■ -Pf^. A manifestly integral algorithm for 
computing any (p^ is provided in [25]. Because (pn has degree <f{n) = HilPi* ~pT~^)i we 
can require each of our polynomials p{x) to be of degree less than ^{n). 

But there is an alternate approach, which is perhaps a little simpler. We represent 
each number by an integer polynomial p{x), and always reduce any exponents modulo 
n (i.e. identify a;" and 1, x'^'^^ and x, etc). When n is even, we can if we like reduce 
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exponents modulo n/2, by identifying a;"'/^ with —1, etc. We equate two polynomials 
p{x), q{x) e Z[x] if 

|p(exp[27riA;/n]) — g(exp[27riA;/n])| < 0.5 

holds for all A;, 1 < A; < n, coprime to n. Of course, '(exp[27rifc/n])'^' is evaluated as 
'cos(27r/cm/n) + i sin(27r/cm/n), to decrease error. The point is that if an algebraic integer 
z (such as z = p(exp[27ri/c/n]) — q'(exp[27riA;/?i])) and all of its Galois associates (j{z) have 
modulus < 1, then z must equal 0. The reason is that ncr(~'^('^)) ^ ^ is the constant 
term in the minimal polynomial of z. 

4-2. Character formulas. The expressions in §3 for S don't depend on the integrality of 
the weight /i, and thus taking the appropriate ratio (4a) gives determinant expressions for 
finite-dimensional Lie group characters. 

4-3. Fusion coefficients. Recall the Kac- Walton formula (2c). Consider any A,// e P^. A 

tensor product coeflBcient T^^ will contribute e(z/) T^^ to the fusion coefficient N^^^ (recall 
that €{iy) and [v] are computed in §3). In particular, any null weight u in the tensor product 
of A and fi can be ignored. When the level is sufficiently high (specifically, Aq + //o ^ k), 
the fusion product of A and n will equal their tensor product. 

4.4- Galois action. By (Ic) we see that the entries of S lie in a cyclotomic field Q[^7v] for 
some root of unity = exp[27ri/A^]. This simply means that each Sxfj_ can be written as 
a polynomial Pxij.{^n) in with rational coefficients. For A^, Bj., . . . , G2, respectively, we 
can take 

AT = 4 (r + 1)k, Ak, Ak, Ak, 12k, 8k, 2k, 2k, 6k 

but usually this is larger than necessary (see [18]). 

Take any Galois automorphism a e Gal(Q[^_;v]/Q)- The group Gal(Q[^_;v]/Q) is iso- 
morphic to the multiplicative group CZ/NZ)^ of numbers coprime to N, so to a we can 
assign an integer i; explicitly, a = ai takes the number Px^{^n) to px^i^^)- example, 
£ = —1 corresponds to complex conjugation. To the Galois automorphism o" is a sign 
CfT '■ P+ ~^ {=t} and a permutation A 1— > cr(A) of P^, such that [16] 

O-('S'am) = 'S'<7(A),M = ^crifJ') Sx,a(n) (9a) 

For example, £ = —1 corresponds to 'charge-conjugation' A t-^ CA and parity e_i(A) = -|-1. 
Some properties and formulas for ecr(A) for A^^^ are given in [20]. For arbitrary a — og, 
the permutation A i— > cr(A) and parity erj{}^ are computed by 

a(A)=[£(A + p)-p] (96) 
e,(A)=e^e(£(A + p)-p) (9c) 

A major motivation for writing this paper was to make this Galois action more ac- 
cessible. In particular, efficient and explicit algorithms for computing \y\ were given in 
§3 for each algebra, and by far the fastest algorithms known to this author for comput- 
ing e{v) use the simplifications to (8) given explicitly in each subsection of §3. The sign 
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is a number-theoretic parity related to quadratic residues, and rarely is needed in the 
applications as it is independent of the weights. It can be computed as follows [16]: 

Write K^l'^s in the form RsfM where i? G Q and M G Z is square- free. Then equals 
the Jacobi symbol (^) G {±1}. This is quickly computed, using standard properties of 
the Jacobi symbol (see e.g. [26]), e.g. quadratic reciprocity and factorisation. 

For example, consider again ^3^5. Here, = 40 works, so for the representatives 
i G (Z/40Z)^ we can take i = ±1, ±3, ±7, ±9, ±11, ±17, ±19. Then, respectively, 

(7,(0) - 0, (2, 2, 2), (2, 0, 2), (0, 6, 0), (0, 6, 0), (2, 2, 2), (0, 0, 0) 
e(^p-p)= 1,1, 1,1, 1,1,1 

€^=1,1,-1,1,-1,-1,-1 

4.5. Modular invariants. The 1-loop partition function of WZW CFT looks like 

z= Yl Mx^^xxx; 

where (among other things) each Mx^ G Z> and MS = SM. Hitting MS = SM with a 
Galois automorphism a implies from (9a) that 

MaX,a,.=Mx^ (10a) 

Ma^^O 6,(A) = €,(//) (106) 

valid for all a. The 'parity rule' (10b) turns out to be especially powerful. Efficient 
algorithms to compute S entries and parities are vital in this context and was the primary 
motivation for writing this paper. See [6] for more details and examples. 

4-6. Branching rules. When X^, C Xg, and the central charge c of Xr,k equals that of 
then we say that we have a conformal embedding. In this case, the level £ modules L{X) 
of xi^^ can be decomposed into a finite direct sum ©jL(A(*'') of level k modules of Xr^\ 
These decompositions are called branching rules. 

All conformal embeddings are known [27] , and most of their branching rules are known 
(see e.g. [28] and references therein), but some branching rules don't seem to appear 
explicitly in the literature, or only appear conjecturally. 

Branching rules can be read off from the appropriate modular invariants (see e.g. 
[29]). In [6] we use the previous algorithms and formulae to obtain branching rules for the 
conformal embeddings. Mostly this merely provides an independent check; however it also 
fills some gaps in the literature. 
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Table 1. Coset representatives for W{D5)\W{E6) and W{AiDq)\W{E7) 



length 


Ci for Eq 


Q for E'j 













Cl = 


/ 


1 




Co — T^ 






C2 = 


ri 


2 










C3 = C2r2 


3 




Ca — C'xT'X 






C4 = C3r3 


4 








C5 = C4r4 , C6 


= C4r7 


5 










= c^rr 


6 








cg = cyre 


, cio = cyrr , 


Cll = C8r3 


7 








C12 = cgrr , 


ci3 = ciors , 


Cl4 = Ciir2 


8 




'-^11' Z 5 '--14 '-'ll' 4 ; 




= Ci2r3 , C16 = Ci3r2 , Ci7 = Ci3r4 , C18 = Ci4ri 
















9 




^lO' i 7 '-'i / ^lo' 4 




= Ci5r2 , C20 = Ci5r4 , C21 = CiqVi , C22 = Ci6r4 


10 




'-'iD'4 7 '-'ly '-^lY'O 




C23 = Cigri , 


C24 = Cigr^ , 


C25 = C2or5 , 










C26 


= C2ir4 , C27 


= C22rs 


11 




^lo' o 7 "-^iy' o 




C28 = C23''4 , 


C29 = C24''3 , 


C30 = C24?'5 , 










C31 


= C26^3 , C32 


= £27^7 


12 




*--20' 2 7 *-"z3 — ^20' D 




C33 = 028^3 , 


C34 = £28^5 , 


C35 = C29r5 , 










C36 = C29r7 


C37 = C3ir2 


C38 = C3ir7 
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C24 = C22r'6 




C39 = C33''2 , 


C40 = C33?'5 , 


C4I = C33r7 , 










C42 = C35r4 


C43 = C35?'7 


C44 = C37r7 
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C25 = ^24^3 




C45 = C39r5 , 


C46 = C39r7 , 


C47 = C40?'4 , 










C48 = C4or7 


C49 = C42r7 


C50 = C44r3 


15 




C26 = C25?^4 




C5I = C45?'4 , 


C52 = C45''7 , 


C53 = C46''3 , 










C54 = C47?'7 


C55 = C49''3 


C56 = C5or4 


16 




C27 = C26?^5 


C57 = 


= C5ir3 , C58 = C5ir7 , C59 = C52r3 , C60 = C53r4 , 










C61 = C54r3 


C62 = C55r2 


C63 = C56r5 



The coset representatives cj are given recursively in terms of the simple reflections r^. 
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Table 2. Coset representatives for W{Ds)\W{Es) 



length 


Ci for Es 









1 




C2 = r-7 


2 






3 




C4 = Cars 


4 




Tc; = r/lT/l = CaTh 


5 




C7 = C5r3 , cg = csrg 


6 






7 

8 




C12 = cgri , ci3 = cgrs , ci4 = Ciors , ci5 = Cure 

'-'Id ^IZ' ; '-'If '-'lo' 5 '-'lo '-'14' 4 5 '-'ly '-'14' 5 '-'ZU '-'lo' / 


9 




C21 = ciers , C22 = Ci7r4 , C23 = cirre , C24 = cigre , C25 = cigry 


10 




C26 = C2ir4 , C27 = C2ir6 , C28 = C22?'3 , C29 = 022^6 , 
^oU ^zo' ( ■) '-'ol '-'24' 5 7 '-'z4' 7 


11 




'-'oo ^2o' 6 5 '-^o4 ^2o' o 5 "-^oo '-'27' 7 5 '-'^ou '-'^zo' b i 
C37 = C29r5 , C38 = C29r7 , C39 = C3ir7 , C40 = Csirg 


12 




— driTo Cao — C'ioTpi Cao — CoATt^ Caa — Co.aT'? Ca^ = Co^aT^ 
^41 ^00' z 7 ^42 ^00' D 7 ^4o '-'o4' O 7 '-'44 ^o4' ( 7 ^4o ^oD' O 7 

'-'4D '-'oD ' ( 7 '-'4 / '-'o / ' ( 7 '-'45 '-'0 / ' 7 '-'49 '-'oy ' D 7 '-'OU '-'oy ' 


13 


C51 


= C4ir6 , C52 = C42r5 , C53 = C42r7 , C54 = C43r7 , C55 = C43r8 , C56 = C45r4 , 

'-'O / ^4:0 ' / 5 '-'OO '-'4i3 ' 5 '-'iDy '-'4 / ' b 5 '-'bU '-'4 r ' c5 5 '-'bl '-'4y ' 


14 




'-^b2 '-'51 ' 7 ^bo '-'51' 7 7 ^d4 '-'52' 4 7 ^b5 '-'52' 7 7 '-'bb '-'52' 7 
^01 *^54' b 7 *^bo '-•54' f? 7 '^'•by '-'5b' 7 7 ^7U '^5b' 7 '^'•71 '-"57' b 7 
C72 = C57r8 , C73 = C59r8 , C74 = C6ir5 


15 


^70 


^OZ' 4 5 ^7o ^DZ' 7 5 ^77 ^DZ' ; ^7o ^d4' 7 5 ^79 ^d4' o 5 ^oU ^DO ' D i 




Cat 


^DO' o ) ^oz ^b7' o ) ^Oo ^d9' O ) ^84 ^d9' o ) ^00 ^7U' ) ^OD Wi' ) 

'--'O ( '--'(o' 5 '--'OO '--'14:' 4 


16 




<^75'd 5 "^9U <--75' 7 5 "^91 <^75' 8 i <^9z <^7d' D 5 "^9oi <^7d'o 5 "--94 <--78'o 7 




'-'9b 


— ''To' 8 ; '^96 — '-'19' 5 ; '-'91 — '^80' 8 ; '^98 — '^82' 5 ; '^99 — '^83' 5 ; '^100 — '^83' 8 ; 
ClOl = C84?'5 5 C102 = 085^6 , C103 = £86^5 , C104 = 087^4 , C105 = £88^3 


17 


C106 = 


= C89r7 , C107 = £89^8 , C108 = Cgore , C109 = C9or8 , Clio = C9ir5 , Cm = C92r8 , 




C112 


= C94r5 , C113 = C94r8 , Cii4 = Cq^T^ , C115 = Cgere , C116 = C97r5 , C117 = C98r4 


18 




C118 = cioare , cng = cioerg , C120 = cio7?'5 , C121 = Ci08r5 , 
C122 = cio8r8 , C123 = Clogrs , C124 = choTq , C125 = cmrs 


19 




Cl26 = Cii8r5 , C127 = Cii8r8 , Ci28 = Ciigr5 , Ci2g = Ci2or4 , C130 = Ci2or6 


20 




Cl31 = Ci26r'4 , C132 = Ci28r-4 , C133 = Ci2gr6 


21 




C134 = Ci33r5 


22 




C135 = Ci34r8 



The coset representatives Cj are given recursively in terms of the simple reflections r^. 



20 



